function dy = iris_EOM(t,y,flag,u)

J = 0.029125;
% J = 0.05;
B = 0.4;
tau = 0.0125;

dy=zeros(2,1);

dy(1)=y(2);
dy(2)=-(1/tau)*y(2)+(B/(J*tau))*u;
end

